A temporally second-order accurate Godunov-type scheme for solving the extended Boussinesq equations     

J.B. Shiach  C.G. Mingham 《Coastal Engineering》2009

A numerical scheme for solving the class of extended Boussinesq equations is presented. Unlike previous schemes, where the governing equations are integrated through time using a fourth-order method, a second-order Godunov-type scheme is used thus saving storage and computational resources. The spatial derivatives are discretised using a combination of finite-volume and finite-difference methods. A fourth-order MUSCL reconstruction technique is used to compute the values at the cell interfaces for use in the local Riemann problems, whilst the bed source and dispersion terms are discretised using centred finite-differences of up to fourth-order accuracy. Numerical results show that the class of extended Boussinesq equations can be accurately solved without the need for a fourth-order time discretisation, thus improving the computational speed of Boussinesq-type numerical models. The numerical scheme has been applied to model a number of standard test cases for the extended Boussinesq equations and comparisons made to physical wave flume experiments.  相似文献   

Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method     

Chao‐Shi Chen  Chien‐Chung Ke  Chia‐Huei Tu 《国际地质力学数值与分析法杂志》2008,32(11):1341-1364

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Resonantly Damped Surface and Body MHD Waves in a Solar Coronal Slab with Oblique Propagation     

I. Arregui  J. Terradas  R. Oliver  J. L. Ballester 《Solar physics》2007,246(1):213-230

The theory of magnetohydrodynamic (MHD) waves in solar coronal slabs in a zero-β configuration and for parallel propagation of waves does not allow the existence of surface waves. When oblique propagation of perturbations is considered, both surface and body waves are able to propagate. When the perpendicular wavenumber is larger than a certain value, the body kink mode becomes a surface wave. In addition, a sausage surface mode is found below the internal cutoff frequency. When nonuniformity in the equilibrium is included, surface and body modes are damped by resonant absorption. In this paper, first, a normal-mode analysis is performed and the period, the damping rate, and the spatial structure of the eigenfunctions are obtained. Then, the time-dependent problem is solved, and the conditions under which one or the other type of mode is excited are investigated.  相似文献   

1998年中国区域降水低频变化的传播机制的初步分析   总被引：7，自引：0，他引：7  

徐国强  朱乾根  薛纪善  何金海 《大气科学》2004,28(5):736-746

使用美国NCEP/NCAR再分析和中国实际降水的逐日资料,研究了1998年夏季中国区域低频降水传播变化的机制.结果显示,(1)中国中东部地区低频降水经向传播都与南海季风和印度季风的变化相对应,低频降水在经向传播中存在不连续现象,是由于低频下沉气流造成的;(2)中国中部区域低频降水纬向传播特征是在不同时间低频降水的纬向传播方向有明显差别,初步认为造成这种现象的可能机制是纬向风气流的变化;(3)高原地区与平原地区之间的低频雨带活动具有不连续性,这是由于在高原东坡存在低频下沉补偿气流造成的.  相似文献   

Observational Zonal Mean Flow Anomalies: Vacillation or Poleward Propagation?          下载免费PDF全文

SONG Jie 《大气和海洋科学快报》2013,6(1):1-7

North-south displacements and meridional vacillations of the eddy-driven jet are widely accepted as the dominant cause of variability of the observational zonal-mean zonal wind anomalies(denoted [u]’).In this study,a new idea regarding the primary variability of the observational [u]’ in the middle latitude troposphere is presented.It is hypothesized that there are two different classes of primary variability of the observational [u]’:the poleward propagation of the [u]’(abbreviated as PP) and meridional vacillations.To validate this hypothesis,one-point correlation maps of [u]’ at 200-hPa during the boreal cold season(November-April) of every year from 1957-2002 are used as a criterion.Twelve PP years,in which the PP events are dominant in the variability of [u]’,and 15 no_PP years,in which the PP events are recessive and the meridional vacillations are dominant in the variability of [u]’,are examined.The results show that the variabilities of [u]’ are different in the chosen PP and no_PP years.In the PP years,the PP events dominate the variability of [u]’;however,the meridional vacillations are prevalent in the no_PP years.  相似文献   

The Relationship Between Abnormal Meiyu and Medium-Term Scale Wave Perturbation Energy Propagation Along the East Asian Subtropical Westerly Jet     

金荣花  杨宁  孙晓晴 《热带气象学报(英文版)》2020,26(2):125-136

The East Asian subtropical westerly jet(EASWJ) is one of the most important factors modulating the Meiyu rainfall in the Yangtze-Huaihe River Basin, China. This article analyzed periods of the medium-term EASWJ variation,wave packet distribution and energy propagation of Rossby waves along the EASWJ during Meiyu season, and investigated their possible influence on abnormal Meiyu rain. The results showed that during the medium-term scale atmospheric dynamic process, the evolution of the EASWJ in Meiyu season was mainly characterized by the changes of3-8 d synoptic-scale and 10-15 d low-frequency Rossby waves. The strong perturbation wave packet and energy propagation of the 3-8 d synoptic-scale and 10-15 d low-frequency Rossby waves are mostly concentrated in the East Asian region of 90°-150°E, where the two wave trains of perturbation wave packets and wave-activity flux divergence coexist in zonal and meridional directions, and converge on the EASWJ. Besides, the wave trains of perturbation wave packet and wave-activity flux divergence in wet Meiyu years are more systematically westward than those in dry Meiyu years, and they are shown in the inverse phases between each other. In wet(dry) Meiyu year, the perturbation wave packet high-value area of the 10-15 d low-frequency variability is located between the Aral Sea and the Lake Balkhash(in the northeastern part of China), while over eastern China the wave-activity flux is convergent and strong(divergent and weak), and the high-level jets are strong and southward(weak and northward). Because of the coupling of high and low level atmosphere and high-level strong(weak) divergence on the south side of the jet over the Yangtze-Huaihe River Basin, the low-level southwest wind and vertically ascending motion are strengthened(weakened), which is(is not)conducive to precipitation increase in the Yangtze-Huaihe River Basin. These findings would help to better understand the impact mechanisms of the EASWJ activities on abnormal Meiyu from the perspective of medium-term scale Rossby wave energy propagation.  相似文献   

SIMULATION OF BOUNDARY LAYER EFFECTS ON A HEAVY RAINFALL EVENT CAUSED BY A MESOSCALE CONVECTIVE SYSTEM OVER THE YELLOW RIVER MIDSTREAM AREA     

赵桂香  王一颉  王晓丽 《热带气象学报(英文版)》2020,26(1):111-124

A heavy rainfall event caused by a mesoscale convective system (MCS), which occurred over the Yellow River midstream area during 7–9 July 2016, was analyzed using observational, high-resolution satellite, NCEP/NCAR reanalysis, and numerical simulation data. This heavy rainfall event was caused by one mesoscale convective complex (MCC) and five MCSs successively. The MCC rainstorm occurred when southwesterly winds strengthened into a jet. The MCS rainstorms occurred when low-level wind fields weakened, but their easterly components in the lower and boundary layers increased continuously. Numerical analysis revealed that there were obvious differences between the MCC and MCS rainstorms, including their three-dimensional airflow structure, disturbances in wind fields and vapor distributions, and characteristics of energy conversion and propagation. Formation of the MCC was related to southerly conveyed water vapor and energy to the north, with obvious water vapor exchange between the free atmosphere and the boundary layer. Continuous regeneration and development of the MCSs mainly relied on maintenance of an upward extension of a positive water vapor disturbance. The MCC rainstorm was triggered by large range of convergent ascending motion caused by a southerly jet, and easterly disturbance within the boundary layer. While a southerly fluctuation and easterly disturbance in the boundary layer were important triggers of the MCS rainstorms. Maintenance and development of the MCC and MCSs were linked to secondary circulation, resulting from convergence of Ekman non-equilibrium flow in the boundary layer. Both intensity and motion of the convergence centers in MCC and MCS cases were different. Clearly, sub-synoptic scale systems in the middle troposphere played a leading role in determining precipitation distribution during this event. Although mesoscale systems triggered by the sub-synoptic scale system induced the heavy rainfall, small-scale disturbances within the boundary layer determined its intensity and location.  相似文献   

Applications of wavelet analysis in differential propagation phase shift data de-noising   总被引：3，自引：0，他引：3  

HU Zhiqun  LIU Liping 《大气科学进展》2014,31(4):825-835

Using numerical simulation data of the forward differential propagation shift （ΦDP） of polarimetric radar,the principle and performing steps of noise reduction by wavelet analysis are introduced in detail.Profiting from the multiscale analysis,various types of noises can be identified according to their characteristics in different scales,and suppressed in different resolutions by a penalty threshold strategy through which a fixed threshold value is applied,a default threshold strategy through which the threshold value is determined by the noise intensity,or a ΦDP penalty threshold strategy through which a special value is designed for ΦDP de-noising.Then,a hard-or soft-threshold function,depending on the de-noising purpose,is selected to reconstruct the signal.Combining the three noise suppression strategies and the two signal reconstruction functions,and without loss of generality,two schemes are presented to verify the de-noising effect by dbN wavelets：（1） the penalty threshold strategy with the soft threshold function scheme （PSS）; （2） the ΦDP penalty threshold strategy with the soft threshold function scheme （PPSS）.Furthermore,the wavelet de-noising is compared with the mean,median,Kalman,and finite impulse response （FIR） methods with simulation data and two actual cases.The results suggest that both of the two schemes perform well,especially when ΦDP data are simultaneously polluted by various scales and types of noises.A slight difference is that the PSS method can retain more detail,and the PPSS can smooth the signal more successfully.  相似文献   

������С����ƽ�����������     

王乐洋  鲁铁定 《大地测量与地球动力学》2014,34(2):55-59

???t??????????ó?????????С????????????????????????????????1????????????????????????2????????????????????????????????????????????Э????????????????1??Э??????????????????????????????????????????????????????????????????????????????????????????????????????????????Ч???  相似文献   

双偏振雷达差分传播相移的五种滤波方法对比分析     

魏庆  ;胡志群  ;刘黎平 《成都信息工程学院学报》2014,(6):596-602

为了对双线偏振雷达的差分传播相移ΦDP进行滤波降噪处理,提高双线偏振雷达识别粒子、估测降水的能力。通过对C波段双线偏振雷达差分传播相移分别采用滑动平均、中值滤波、FIR滤波、卡尔曼滤波和小波分析5种不同方法地介绍,利用中国气象科学研究院灾害天气国家重点实验室移动式C波段双线偏振雷达(POLC)2013年外场试验观测资料,对比分析了经5种方法处理后的差分传播相移ΦDP径向数据,及其分别在降水估测中应用的效果。结果表明:采用FIR滤波方法、卡尔曼滤波方法和小波分析方法对ΦDP径向数据进行处理的效果比传统的滑动平均和中值滤波方法好;其中小波分析处理的ΦDP径向数据波动最小,再经过变距离法拟合的KDP数据的负值能够得到最有效地抑制,最后应用到R(KDP)降水估测中的精度最高。  相似文献